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ABSTRACT 

The synchrotron self-Compton (SSC) emission from Gamma-ray Burst (GRB) forward shock can 
extend to the very-high-energy (VHE; E-y >100 GeV) range. Such high energy photons are rare and 
are attenuated by the cosmic infrared background before reaching us. In this work, we discuss the 
prospect to detect these VHE photons using the current ground-based Cherenkov detectors. Our 
calculated results are consistent with the upper limits obtained with several Cherenkov detectors 
for GRB 030329, GRB 050509B, and GRB 060505 during the afterglow phase. For 5 bursts in our 
nearby GRB sample (except for GRB 030329), current ground-based Cherenkov detectors would not 
be expected to detect the modeled VHE signal. Only for those very bright and nearby bursts like 
GRB 030329, detection of VHE photons is possible under favorable observing conditions and a delayed 
observation time of <10 hours. 
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1. INTRODUCTION 

On February 28, 1997, the first X-ray afterglow of 
a Gamma-ray burst (GRB) was detected, leading to 
the ide ntification of its p rogenitor at cosmological dis- 
tances (jCosta et al.l 119971 ). In a few days, the afterglow 
faded away with time as a power law. This behavior is 
satisfactorily explained in the spherical (isotropic) fire- 
ball model involvin g relativistic ejecta decel erated by cir- 
cumburst medium (Mcszaros & Rccs 1997). The intro- 
duction of collimated jets relaxes the energy requirement 
of GRBs by a factor of several hundred, as we ll as ex- 
plains the steeper tem poral decay of afterglows (jRhoadsl 
|]^9lt[SarTet^fT999l) . 

GRBs are extra-galactic sources of GeV and proba- 
bly higher energy photons. Evidences of distinct high- 
energy component have been accumulated by EGRET 
onboard the Comp ton Gamma- Ray Observatory: (1) 
iHurlev et al.l (|1994[ ) reported the detection of long- 
duration MeV-GeV emission of GRB 940217 lasting up 
to 1.5 hour after the keV burst, including an ~18 GeV 
photon. This burst is the longest and the most energetic 
among th ose GRBs with detecte d high-energy emission 
so far; (2) iGonzalez et a l. (2003) revealed a high-energy 
component of GRB 941017 temporally and spectrally dif- 
ferent from the low-energy component. 

In the fireball model, synchrotron emission of shock- 
accelerated electrons is commonly thought to produce 
prompt 7-ray emission as well as afte rglow emission 
at lower energies (e.g., ISari et all [19 98). It is natural 
to expect that these photons are inverse-Compton up- 
scattered by electrons, giving rise to a higher e nergy com- 
ponent peaking at sub-GeV to TeV energies (jWei fc Lul 
ll998tlSari fc Esinll200l[) . When electrons scatter the self- 
emitting synchrotron photons, synchrotron self-Compton 



(SSC) emission is resulted. In the external shock sce- 
nario, the temporal profile of the SSC emission from for- 
ward shock electrons is similar to that of the low en- 
ergy afterglow emission and no significant time lag is 
expected. 

The Swift satellite, thanks to its rapid response time 
and accurate localization, has started a new era of re- 
search on GRBs. Different modifications to the standard 
afterglow model are put forward to explain the pecu- 
liar behaviors exhibited in the X-ray l ight curves, in par- 
ticular the shallow declining phase (|Zhang et all 120061 : 
iNousek et aLll200fih . Recently, the SSC emission of the 
modified forwa rd shock has been extensively discussed in 
the literature (IWei fc Fanl 120071: iGou fc Meszarod 120071 : 
iFan et alJfeOOa I G alii fc Pird 120071 : lYu et all 120071) a nd 
applied to the case of GRB 940217 (|Wei fc Fanll2007f ). 

The AGILE (Astro-rivelatore Gamma a Immagini 
L'Eggero) satellite, launched on April 23, 2007, is ded- 
icated to high-energy 7-ray astronomy. The Fermi 
Gamma-ray Space Telescope (FGST) was launched on 
June 11, 2008. The Large Area Telescope (LAT) on- 
board covers the energy range from 20 MeV to 300 
GeV and its effective area is about 5 times larger than 
that of EGRET at GeV energies. The first GRB ob- 
servations with LAT have resulted in detection of pho- 
tons w ith energies larger than ~ 1 GeV from several 
GRBs JOmodeill2008t iTaiima et al.ll2008fl . iDermer et al 



(2000) . IZhang fe Meszarosl (|2001bf ). and IWang et al 



20011 ) predicted promising and detectable SSC emission 
from the forward shock with FGST out to z~l. 

Most of the discussions in the literature have focused 
on the afterglow emission from tens of MeV to GeV. 
LAT can also detect very-high-energy (VHE; >100 GeV) 
afterglow emission. However, with a small effective 
area ~ 10 cm 2 , it is very hard to have a significant 
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Cherenkov telescopes such as H.E.S.S. 1 , MAGIC 2 , and 
VERITAS 3 may serve better at energies above ^100 GeV 
because of their much larger effective area (~ 10 s — 
10 9 cm 2 ) and a high rejection rate of hadronic back- 
ground. The effective collectin g area of Cherenkov tele - 
scopes increases with energy jAharonian et al.1 l2006bh . 
Some of these large area Cherenkov detectors have been 
used to set constraints on the poss i ble VHE afterglow 
component of GRBs (I Albert et all 120071 : iHoran et all 
120071: lAharonian et alj|2009t ). It is thus desirable to see 
whether these results are consistent with the predictions 
of the fireball model. Our aim of this paper is also to 
investigate the prospect of significant detections in the 
future. To have a reliable estimate of the afterglow emis- 
sion at energies above 100 GeV, one needs to calculate 
the forward shock emission (both synchrotron and SSC 
emission of the shocked electrons) carefully. The attenu- 
ation of VHE photons by the cosmic infrared background 
is also taken into account. Since the attenuation effect for 
photons with an energy >100 GeV is more severe for 
high-redshift GRBs, we limit our GRB sample to nearby 
events. 

This paper is organized as follows: in Sj2j we describe 
the GRB afterglow model, introduce the code that is 
used in the afterglow modeling and the calculation of 
the SSC emission from GRB forward shock. In we 
present the expected results of the SSC model using rea- 
sonable parameter values for GRBs. In £01 we describe 
the GRB sample which includes six nearby GRBs with 
sufficient multi-wavelength afterglow data and predict 
their corrected energy flux after the attenuation by the 
cosmic background during the afterglow phase, which is 
then compared with the available observational data. We 
summarize our results and discuss their implications in 
SJS] We conclude in [JS] 

2. AFTERGLOW MODELING 
2.1. GRB Afterglow Model 

While synchrotron emission is widely considered to 
be responsibl e for the radio, optical, and X-ray af- 
terglows (e.g. iSari et aLI 119981 ). inverse Compton scat- 
tering of forward shock phot ons i s considered in de - 
tails bv lWei fc Lul (fl998l l2000h and lSari fc E"sml (|200ll) . 
Inverse Compton scattering may considerably change 
the temporal and spectral behavior of GRB afterglows, 
and its cooling effect on electrons accelerated in ex- 
ternal shocks will contribute to th e photon spectra 
at sub-GeV to TeV ener g y range jM eszaros fc Reed 
1994 iDermer et al.1 120001: IZhang fc Meszarosl l2001bl : 
Wang et al.ll2001ri . 

In the afterglow model, both synchrotron emission and 
inverse Compton emission are taken into account. It is 
assumed that: (1) the external medium is homogenous 
with a density n or a wind profile n oc R~ 2 ] (2) the rela- 
tivistic jet is uniform, i.e. energy per solid angle is inde- 
pendent of direction within the jet; (3) the shock parame- 
ters (e e and e b , fractions of the shock energy given to the 
electrons and the magnetic field, respectively) are con- 
stant; (4) the energy distribution of electrons accelerated 
in shocks follows dN e /dE oc E~ p ; (5) the possible achro- 

1 http://www.mpi-hd.mpg.de/hfm/H.E.S.S./H.E.S.S..html 

2 http://wwwmagic.mppmu.mpg.de/ 

3 http:/ /veritas. sao.arizona.edu/ 



matic flattening in the afterglow light curve is due to 
energy injection in the form Ek oc t 1 ~ q (|Cohen fc PiranI 
fl999ElZhang fc Meszarodl20^ll or E k oc [1 + (t/T) 2 ]~ 1 
with T being the initial spin-down time scale (|Dai fc Lul 
fl99l . 

The parameters involved in this afterglow model in- 
clude: E (the initial isotropic outflow energy), O (the 
initial half-angle of the jet), n (the density of the homo- 
geneous external medium) or A (the wind parameter), 
p (the power-law index of energy distribution of shock- 
accelerated electrons), e c , and ee (shock parameters). In 
the case where energy injection is necessary, three addi- 
tional parameters: L e je (the injected luminosity in the 
rest frame), the timescale of energy injection and q, are 
included. 

2.2. A Brief Description of the SSC Model 

The code used in our afterglow modeling and the 
pre diction of the SSC emission is that developed 
by iFan et ail (|2008l ). who carry out numerical calcu- 
lations of synchrotron and SSC emission of the exter- 
nal forward shock in the afterglow phase 4 . The re- 
verse shock emission, predicted i n the fireball model 
but not detected i n most events (jRoming et al. I [20061 : 
iRvkoff et al. Il2009l ). is not take n into account. 

The key treatments (see §3 of IFan et al.ll2008l . for de- 
tails) are as follows: (i) The dynamical ev olution of the 
outflo w is followed using the formulae in iHuang et al.l 
(2000), which describes the hydrodynamics in both rel- 
ativistic and non-relativistic phases, (ii) The arbitrary 
assumption that the distribution of shocked electrons is 
always in a quasi-stationary state is considered to be un- 
satisfactory and the energy distribution of electrons is 
calculated by solving the continuity equation with the 
power-law source fu nction Q = K7J P , no rmalized by a 
local injection rate (jModerski et al.l l2000). (iii) The ob- 
served flux is integrated over the "equal-arrival surface" . 
(iv) The Klein-Nishina correction of the inverse Comp- 
ton emission has been included, (v) Energy injection into 
the outflow is considered necessary in reproducing some 
multi- frequency afterglow data of some GRBs. This may 
change the dynamics significantly. 

By fitting the low-energy multi- waveband afterglows, 
from radio to X-ray band, parameters involved in the 
afterglow model are gotten. Simultaneous afterglow data 
in at least two well- separated wavebands are needed to get 
a relatively well- constrained set of parameters. 

A rough estimate of the energy-integrated VHE after- 
glow flux (without correction by the cosmic background) 
is given by 

F>iooGeV oc (1 + *l Lssc max{(0^, Kf)^} (1) 

L 

where L ssc is the tot al luminos i ty of the SSC emis- 
sion (see eq.(23-27) in IFan et all 120081 for the expres- 
sion), v^f and are the typical SSC emission fre- 
quency and the SSC cooling frequency of the forward 

4 The cooling of the forward shock electrons by the so-called 
"central engine afterglow" emission (the emission powered by the 
re-activity of the central engi ne in the afterglow phase, like the 
flares (e.g., Nousek et al. 2006) or the plateaus followed by a sud- 
den drop (e.g., Zhang 2009)) and the corresponding inverse Comp- 
ton emission can also be calculated self-consistently. But for our 
current sample people did not see such emission. 
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Fig. 1. — Temporal evolution of the observed high energy spec- 
trum of SSC afterglows. The three dotted lines are the spectra at 
different epochs. All spectra, starting from (top) 0.5 hour, 2 hours, 
and (bottom) 10 hours, are integrated over 0.5 hour. They are cal- 
culated with the following parameter values: Eq =5 X 10 51 erg, 
O = 0.4, n = 1.0cm -3 , p = 2.2, e c = 0.3, e B = 0.01 and 
z = 0.16. The solid and dashed line are calculated with the 
same parameters represented above during the first epoch but 
occurring at larger redshift, z = .5 and z = 1.0, respe ctively. 
The sensitiv ity curves of EGRET HThompson et al.|[T993fi, AG- 
ILE/ GRID (G alli fc Guctta 200|), FGST/LAT l IGalli fc Guettal 



[20081. VERITA S l|Horan et al.l 120071 ) and HESS (assuming F 
2.6, Tarn 2009) for an integration time of 0.5 hour are plotted as 
labeled. Magic 2-sigma upper limits derived from 30-min obser- 
vations of GRB 060206 are also plotted, taken from Albert et al. 
(2007). 



shock electrons (see eq. (33-34) in iFan et al.l 120081 the 
case of k = 0, for the expressions) and Dl is the lu- 
minosity distance of the event. 

3. MODEL PREDICTION 

High energy photons, especially those in the TeV 
range, will be attenuated by the cosmic background 
light. Various models of the spectral energy dis- 
tributio n of the cosm i c infrared background are pro- 
posed (iPrimack et al] I200H PTbtani fc Takeuchl 120021 : 
iKneiske et"aH 120021 : fStecker et all 120061 ). but all these 
models give comparable opacities for low redshifts. In 
this work, a level consistent with a study of two 
distant blazars and galaxy counts is used (P0.45, 
lAharonian et al.ll2006aft . 

We adopt reasonable values of parameters for nearby 
GRBs and predict the spectra in high-energy to VHE 
range. After corrected for the attenuation by extragalac- 
tic background, we compare them with the sensitivity 
levels of 7-ray instruments. 

Parameters assumed and the time-averaged spectra, 
including both synchrotron and SSC components from 
the forward shocks, are shown in Figure [TJ For this ficti- 
tious burst, current Imaging atmospheric Cherenkov tele- 
scopes would be more likely than space satellites such as 
EGRET, AGILE/GRID and FGST/LAT to detect the 
modeled emission, as seen in Figure [T] 



4. VERY HIGH ENERGY AFTERGLOW EMISSION OF 
NEARBY GRBS 

For photons with energy higher than ~ 100 GeV, the 
attenuation due to interaction with background photons 
is significant if the source has a high redshift. There- 
fore nearby bursts (those with z<0.25) are chosen in this 
study. 

4.1. The GRB Sample 

To predict the VHE afterglow emission of nearby GRBs 
and compare the calculated results with sensitivity of 
different 7-ray telescopes, GRBs in our sample must sat- 
isfy: nearby events to alleviate the attenuation effect; at 
least two independent waveband low-energy afterglows 
are recorded to get relatively constrained parameters; the 
low-energy afterglow can be reproduced according to the 
afterglow model. 

In this work, we consider nearby GRBs (z < 0.25) with 
relatively high luminosity and multi-wavelength after- 
glow data sufficient to meaningfully constrain the prop- 
erties of the GRBs (i.e. the model parameter values as 
described in fl2~Tj) up to March 2007. Five GRBs meet 
such criteria: GRB 030329, GRB 050509B, GRB 050709, 
GRB 060505, and GRB 060614. Though having a rel- 
atively large redshift of z~0.55, GRB 051221A is also 
considered in this work because it is one of the brightest 
short GRBs detected so far. 

GRB 030329 tri ggered the High Energ y Transient 
Explorer, HETE-2 (|Vanderspek et al] 120041 ) . Very de- 
tailed BVRI afterglow light curves, spanning from ~0.0 5 
to ~ 80 day s , wer e compiled by iLipkin et al.l |2004). 
iTiengo et all (|2004ft reported XMM-Newton and Rossi- 
XTE late-time observations of this burst. Based on the 
emission and absorption lines in the optic al afterglow, a 
redsh ift of z=0.1685 has been identified (|Greiner et al] 
l2003h . 

The X-Ray Telescope (XRT) onboard Swift began ob- 
servations of GRB 050509B 6 2s after the trigger of the 
Burst Alert telescope (BAT) ([Gehrels et al 112005ft. Op- 
tical a nd infrared data were repor ted in Bloo m et al] 
(|2006ft . iProchaska et all ([20051 ) and lBloom et all (|2005l ) 
reported a redshift of z~0.22 based on numerous absorp- 
tion features and a putative host galaxy, respectively. 

GRB 050709 was discovered by HETE- 
2 (jVillasenor et al.l 120051 ). Its prompt emission lasted 
70 ms in the 3-400 keV energy band, followed by a 
weaker, soft bump of ^100 s duration. The optical 
counterpart of this burst was observed with the Danish 
1.5-m telescope at the La Silla Observatory. The 
observations started 33 hours af ter the burst an d 
spanned over the following 18 days (jHiorth et al.ll2005T ). 
Observations with the Chandra X-ray observatory 
revealed a faint, uncata logued X-ray s ource inside 
the HETE-2 error circle (|Fox et al] [2005ft . which was 
coincide nt with a pointlike object embedded in a bright 
galax y (|Jensen et al] [2005) at z = 0.16 (|Price et al] 
[200l . 

GR B 051221A was localized by BAT (jCummings et al.l 
2005) and also simultaneously observed by the Konus- 



Wind instrum ent 
[Burrows et al] [200 



The X-ray (~ 10 2 
and t he optical 



2 x 10 b s; 



10 4 



4 x 



10 5 s; ISodcrbc rg et al.l l2006ft afterglow light curves of 
GRB 051221A were well detected, while in the radio 
band only one detection followed by several upper lim- 
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its are available (|Soderberg et al.|[2006f ) . ISoderberg et al.1 
(2006) detected several bright emission lines, indicating 
a redshift of z = 0.5464. 
GRB 060 505 was detected b y BAT in the 15-15 
((Palmer et all 120061 : iHullinger et"aT1 l2006f ). 
2006h reported the detection of the opti- 
later confirm ed by VLT FORS2 obser- 
vations (jThoene et al.ll2006l) . XRT detected a source 
which wa s located about 4" from a galaxy with 
z=0.0894 (jConciatore et "aT1l2006fh 



keV band 
lOfek et all 
cal tran sient, 



GR B 060614 triggered both Swift-BAT ([Parsons et al.l 
I2006T ) and Konus-Wind dGolenetskii et al.ll200fif ). XRT 
found a very bright (^1300 counts s _1 ) un-catalogued 
source inside the BAT error circle. Ground-based 
optical and infrared follow-up obser vations were per- 
formed using several instruments (e.g.. lCobb et eil.ll2006t 
iSchmidt et al.1 120061 ). Based on the detection of the 
host galaxy emissi on lines, a redshi ft of z = 0.125 
proposed by [Price et al.1 |2006) and confirmed 
Noted that the classification 
the commonly-used 
long/short burst schem e, since it has a long duration but 
no ac companying SN ijGehrels et all [20061 iFvnbo et all 
2006). 



was 



by \fnce 
bv lFugazza et al.1 (|2006l ). 
of GRB 060614 is ambiguous 



4.2. Constraining the Model Parameters 

The available multi-frequency afterglow data are then 
used to obtain the model parameters. In this work, we 
have reproduced the multi-frequency afterglow data of 
GRB 030329 and GRB 060614. 

The well-sampled distinguishing afterglow behavior of 
GRB 030329 has gained much attention. Some authors 
concentrated on the rebrightening occurring at 1.6 days 
after the trigger and considered different mechanisms to 
explain the rebrightening f eatures seen in the optical light 
curves (|Huang et al.ll200H ) . We c oncentrate on the multi- 
wave band emission from radio (jBerger et al.ll2003T). op- 
tical (jLipkin et al.l f2004h to X-ray band (jTiengo et all 
2004) for the purpose in this paper. We show in Fig- 
ure [2] that, with a certain set of parameters, the numer- 
ical results can describe the observed data in all three 
wavebands. Fluctuations were captured in R band af- 
terglow light curves after 5 x 10 4 s from the burst trig- 
ger, which may imply the multiple energy injection into 
the outflow dHuang et al.ll2006h or a two component jet 
(jBerger et al.ll2003[ ). We ignore these details and focus 
on the general trend of the optical emission (particularly, 
in our calculation the energy of the relativistic ejecta is 
a constant). As shown in Figure [H the difference is that 
in the time range 5 x 10 4 — 10 5 sec our approach gives a 
(little bit) brighter optical emission, so will be the high 
energy emission. 

The modeled and observed afterglow light curves of 
GRB 060614 are shown in Figure [3] Energy injection, 
starting around 30 minutes after the GRB onset, is 
needed in the afterglow modeling to reproduce the in- 
crease in flux (instead of simple power-law decay seen in 
other GRBs) . The early X-ray flux before 500s from the 
burst trigger, which is much brighter than the modeled 
flux, results from the dominating contribution from the 
prompt emission. 

Table [1] lists the physical parameters derived from 
the afterglow modeling for GRB 030329 and GRB 
060614. Parameters of GRB 050509B, GRB 050709, 
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Fig . 2.— GRB 03 0329 afterglow data in the 15 GHz l IBerger et al.l 
I2003T), -R-band JLipkin et al.l I2004T) and 0.5-2keV band 
UTiengo et al.l I2004j y! Symbols indicate data points as la- 
belled. Solid lines represent the modeled i?-band and 15 GHz 
emission (left plotting axis) and the dotted line represents the 
modeled X-ray emission (right plotting axis), respectively. 
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Fig. 3.— G RB 060614 afte rglow data in the i?-band and X-ray 
band (see also Xu ct al. 2009). Crosses represent data recorded by 
Swift-XRT (0.2-10 keV), and circles represent i?-band data. Solid 
and dashed lines represent the modeled i?-band and 1 keV emission, 
respectively. 

GRB 051221A and GRB 060505 are taken from the lit- 
erature, also listed in Table 1. 

4.3. VHE Gamma-ray Observational Data 

We are interested in VHE observations during the af- 
terglow phase when the SSC is likely to dominate (see 
iJS]). VHE 7-ray afterglow data of three of the GRBs 
in the sample (i.e. GRB 030329, GRB 050509B, and 
GRB 060505) are available. 
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4.3.1. GRB 030329 

iHoran et al.l (|2007l ) reported a total of 4 hours of obser- 
vations, which spanned five nights, using the Whipple 10- 
m telescope. No evidence for VHE 7-ray signal was found 
during any of the observation periods. When combining 
all data, a flux upper limit of 1.4 x 10 _11 ergcm _2 s _1 was 
derived. The first observation, lasting for about an hour, 
was started 64.6 hours after the burst. The 99.7% c.l. 
flux upper limit above an energy of ~ 400 GeV derived 
from this observation is shown in Table O as well as in 
Figure HJ 

The 28-minute H.E.S.S. observatio n of GRB 03032 9 
were taken 11.5 days after the burst (|Tam et al.ll2008f ). 
Since the burst position was located above the northern 
hemisphere, the zenith angle of the GRB observation was 
relatively large, i.e. 60°, resulting in an energy threshold 
of 1.36 TeV. No evidence for VHE 7-ray signal was found. 
The 99% c.l. flux upper limit (> 1.36 TeV) is 3.4 x 
10 _11 ergcm~ 2 s~ 1 , assuming a photon index of T = 3. 

4.3.2. GRB 050509B 

The observations of this burst using the STACEE de- 
tector employ an " on-off observation mode" and contain 
two 28-minute on/off pairs. The first on-source observa- 
tion started 20 minutes after the burst and the second 
80 minutes after the burst. After data quality cuts, about 
18 minutes of useful on-source data remain in each ob- 
servation. No evidence for VHE 7-ray signa l above the 
energy threshold of 150 GeV was reported bv I Jar vis et al.l 
(2008). The 95% c.l flux upper limits (above 150 GeV, 
assuming a photon spectrum of dN/dE ~ E~ 2 - 5 ) were 
3.8 x 10 _10 ergcm _2 s~ 1 and 4.5 x 10 _10 ergcm _2 s _1 for 
the first and second on-source observation, respectively 
(A. Jarvis private communication). 

4.3.3. GRB 060505 

The H.E.S.S. observations began 19.4 hours after the 
burst and lasted for 2 hours (|Tam et al.ll2008l ). No ev- 
idence for VHE 7-ray signal was found. The 99% c.l. 
flux upper limit (> 0.45 TeV) is 8.8 x 10 _12 ergcm _2 s -1 , 
assuming a photon index of T = 3. 

4.4. Comparison to Observations 

Based on the parameters obtained in §4.2, the GeV- 
TeV emission is obtained using the code described in 
§2.2. 

We depict the calculated high energy afterglow spec- 
trum in Figure [U which shows the time-integrated high 
energy afterglow spectrum of these six events. The solid 
and dashed lines represent the intrinsic SSC spectra and 
corrected spectra for each GRB, respectively. The ab- 
sorption is bas ed on the cosmic infra red background 
model "P0.45" (jAharonian et alJl2006aT > 5 . Such model 
is constrained by the upper limits provided by two un- 
expectedly hard spectra of blazars at optical/NIR wave- 
lengths and is close to the lower limit from integrated 
light of resolved galaxies. 

In order to compare with the VHE observational data 
which are usually given in integrated photon fluxes, we 
integrate the spectra above the energy threshold. We 

5 This implies a gamma ray horizon at a redshift of about 0.2 
(0.05) for 500 GeV (10 TeV) gamma rays. 



consider first the GRBs with VHE data. These include 
GRB 030329, GRB 050509B, and GRB 060505. In Ta- 
ble 2 we list the modeled integrated energy fluxes after 
correction due to interaction with photons from cosmic 
infrared background, as well as the VHE 7-ray observa- 
tions and the derived upper limits. All predicted fluxes 
are below the upper limits derived from the VHE obser- 
vations. 

We then investigate whether a sensitive VHE instru- 
ment is expected to detect the predicted VHE signal from 
nearby GRBs during the late afterglow phase. We use 
H.E.S.S. sensitivity as an example of an array of sen- 
sitive atmospheric Cherenkov telescopes. The sensitiv- 
ity level of H.E.S.S. detec tor is shown in Figure El as- 
suming a r=2.6 spectrum (Aharonian et al. 2006b). As- 
suming softer spectra, the level is higher, and the dif- 
ference is about 50% between T=2.0 and T=3.0 (c.f. 
I Aharonian et al.ll2005[ ). We show the temporal evolution 
of energy fluxes (>200GeV) of six GRBs in our sam- 
ple in Figure [5l indicating only the VHE signal from 
GRB 030329 may be above the H.E.S.S. sensitivity. For 
GRB 030329 which is a bright burst with low redshift, 
the expected energy flux would be high enough to be de- 
tected with a delayed observation time of < 10 hours if 
GRB position was favorable, i.e. with zenith angle < 20° 
(and thus an energy threshold of ~ 200 GeV is attained) . 

5. DISCUSSION 

In this paper, we have calculated the SS C emission 
from the forward shock electrons following iFan et aT] 
(2008). We shall discuss here the importance of other 
radiation processes in the late afterglow phase. 

Possible VHE 7-ra y emiss i on in i tiated from protons 
has b een suggested (jTotanil 119981 : iBottcher &; Dermerl 
1998). However, the proton-synchrotron component, as 
well as the hadron-related photo-meson electromagnetic 
components, is in most cases overshadowed by the SSC 
component of electrons in the afterglow phase. This is 
especially the case for the parameter values of e e and 
€b used here in the mod eling of these six GRBs (e.g., 
IZhang fc Meszarosl[2001bf l. 

Additional IC component will play a role and will en- 
hance any VHE emission. SSC is, in general, only a lower 
limit. Another possible contribution to VHE emission is 
related to the "central engine afte rglow" emission, like 
the flares (e.g.. iFalcone et aLll2007l) or the p lateaus fol- 
lowed by a sudden drop (e.g.. IZhang 1120091 ) . The SSC 
radiation of the late internal shocks or the external in- 
verse Compton radiation in the external forward shock 
front may be able to give rise to some VHE emission sig- 
nals (see Fan & Piran 2008 for a review). However, no 
"central engine afterglow" emission has been reported in 
the late afterglow phase for the six bursts we studied. So 
it is hard for us to estimate its ability to enhance the 
detectability of the VHE emission. 

As shown in Figure EJ we only expect detectable sig- 
nal using a sensitive ground-based 7-ray detector for a 
bright and nearby GRB like GRB 030329. The rate 
of such nearby and energetic GRBs is very uncertain. 
GRB 940217 might h ave been another event of this 
kind (|Wei fc Fanll200l . 

Several factors that reduce the chance of detecting 
VHE photons can be summarized as follows: Firstly, as a 
result of large zenith angles (e.g., 60° for GRB 030329), 
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Fig. 4. — Modeled time-integrated 0.1 GeV — 20 TeV afterglow spectra of six GRBs, in comparison with VHE upper limits (triangles). 
Dotted and solid lines represent the spectra with and without correction by the cosmic infrared background, respectively. For GRB 030329, 
GRB 050509B, and GRB 060505, the spectra were integrated over the corresponding time intervals during which the upper limits were 
derived, as shown in Table [2] For GRB 030329, thick (upper) lines indicate the modeled spectrum for the Whipple observation time, 
and thin (lower) lines for the H.E.S.S. observation time. The data points are plotted at the corresponding average photon energies. The 
modeled spectra of the remaining three bursts are obtained by integrating the spectra over a time period of 2 hours, starting from 10 hours 
after the trigger. 



7 



1(T 



GRB 030329 



JHESS sensitivity *■ - 




10 10 
time (second) 



Fig. 5. — The temporal evolution of modeled VHE integral en- 
ergy fluxes above 200 GeV for GRB 030329 (black dashed line), 
GRB 050509B (blue solid line), GRB 050709 (green solid line), 
GRB 051221A (red dash-dot-doted line), GRB 060505 (blue dot- 
ted line) and GRB 060614 (magenta dashed line). The red solid 
line represents the HESS sensitivity if the observations start right 
at the GRB onset. 

the energy thresholds of some observations are relatively 
high (~1.4 TeV). Any VHE photon is severely attenu- 
ated by the cosmic infrared background, unless the level 
is very low. Secondly, the observations were taken at very 
late epochs, e.g. 11.5 days after the burst for H.E.S.S. 
observations of GRB 030329, when expected VHE flux 
had largely decayed. Thirdly, the fraction of low-redshift 
GRBs is small, e.g. see the catalogue in iButler et al. I 
(|2007f > . For GRB 051221A (at z = 0.55) studied here, 
the attenuation is severe at energies >200 GeV. 

Despite these practical limitations, detection of VHE 
afterglow emission of GRBs is probable. Those GRBs 
close enough (z<0.5) and with an intrinsic high luminos- 
ity (like GRB 030329), can be detected above ~200 GeV 
when the observation is taken within ~10 hours after the 
burst. From Eq. 1, GRBs with low z, large Eq, large e e , 
and small e# are more likely to be detected in the VHE 
band. Besides, in general the expected VHE afterglow 
flux decays as t~ a where a is at least one. Thus, obser- 
vations with shorter delay are more likely to probe the 
predicted VHE emission. 

The main purpose of this work is to provide a rela- 
tively reliable prediction of the detectability of the VHE 
emission from nearby GRBs. It is also interesting to ex- 
plore the role of the VHE afterglow emission detection 
in revealing the GRB physics. Perhaps the most robust 



conclusion is that the VHE emission in late afterglow 
phase cannot be attributed to the synchrotron radiation 
of the forward shock, for wh ich the maximal pho ton en- 
ergy is w 30r/(l + z) MeV (iCheng fc Weilfl99l) . where 
r is the bulk Lorentz factor of the decelerating outflow. 
In general, for es > £e, the inverse Compton emission 
is expected to be very weak. So the detection of VHE 
afterglow emission from GRBs requires that £b "C e e , 
in support of the current radio/optical/X-ray afterglow 
modeling (e.g., see Tab.l). Together with Swift, AG- 
ILE, FGST and some ground-based optical telescopes, 
ground-based 7-ray detectors can provide us continuous 
spectra in the optical to TeV energy band during the af- 
terglow phase. A self-consistent modeling of these data 
in a very wide energy band, in principle, will impose 
very tight constraint on the physical parameters and on 
the environment. However, given the small number of 
the VHE photons expected from a single burst, we do 
not expect that the VHE emission can help us a lot to 
achieve such goals. 

6. CONCLUSIONS 

In this work, we discuss the prospect of detecting VHE 
7-rays with current ground-based detectors in the late 
afterglow phase. During this phase, the dominant radia- 
tion process in the VHE 7-ray regime is the SSC emission 
from the forward shock electrons. Klein-Nishina effects 
and attenuation by the cosmic infrared background, both 
known to suppress the VHE 7-ray spectra, were taken 
into account. To minimize the attenuation effect, we 
chose a sample of six nearby GRBs in this study. We 
have calculated the detailed SSC emissi o n num erically 
using the model developed bv iFan et al.1 (|2008f ). with a 
set of parameters which are able to reproduce the avail- 
able multi-wavelength afterglow light curves. The re- 
sults are consistent with the upper limits obtained using 
VHE observations of GRB 030329, GRB 050509B, and 
GRB 060505. Moreover, assuming observations taken 
10 hours after the burst, the VHE signal predicted from 
five GRBs is below the sensitivity level of a current sen- 
sitive atmospheric Cherenkov detectors mainly due to 
the low fluence of these outflows. For those bright and 
nearby bursts like GRB 030329, a VHE detection is pos- 
sible even with a delayed observation time of ~10 hours. 
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TABLE 2 

VHE GRB OBSERVATIONS AND MODEL PREDICTIONS 



GRB 


telescope 


Tobs - Tqrb a 


exposure 


energy threshold 
(GeV) 


energy flux upper limit set 
by observations (erg cm -2 s — 1 ) 


predicted energy flux 
(erg cm -2 s — x ) 


references 


030329 


H.E.S.S. 


11.5 days 


28 min 


1360 


3.4 X 10- 11 


8.5 X lO" 15 


1 


030329 


Whipple 


64.55 hours 


65.2 min 


400 


5.8 x lO" 11 


6.7 x 10~ 13 


2 


050509B 


STACEE 


20 min/80 min 


28 min 


150 


3.8xl0~ 10 /4.5 x 10~ 10 


2.2xl0 _16 /5.4xl0~ 17 


3 


060505 


H.E.S.S. 


19.4 hours 


2 hours 


450 


8.8xl0~ 12 


2.5xl0~ 15 


1 



References. — UTam et ahl |(2008T) ; 2 IHoran et"aTI J2007h : 3 IJarvis et all ([20081 1 
a The time between the start of the GRB and the beginning of observations for different telescopes. 



